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OBJECTIVE 


Develop  and  analyze  advanced  signal  processing  techniques  that 
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will  provide  improved  performance,  reliability,  and  maintainability 
while  reducing  costs  of  signal  processing  systems  for  radar,  sonar, 
and  electronic  warfare. 

RESULTS 

1.  An  exact  computational  method  is  given  for  determining  the 
detection  threshold  to  be  used  when  a square-law  integrator  is  applied 
to  the  output  from  a nonrecursive  moving  target  indicator  (MTI). 

2.  A table  of  bias  levels  is  given  for  a two-pulse  (single- 
delay) MTI  for  several  values  of  the  number  N of  pulses  integrated 

and  the  probability  PFA  of  a false  alarm. 

I 

3.  The  characteristic  function  of  the  output  from  the  square- 
law  integrator  is  determined  for  the  case  where  the  MTI  input  contains 
a signal . 


RECOMMENDATION 

Continue  the  analysis  of  the  effect  of  MTI  on  incoherent  integration. 
Specifically,  determine  exact  (or  approximate)  methods  for  computing 
detection  probabilities,  and  investigate  optimum  methods  for  constant 
false-alarm  rate  (CFAR)  processing  for  an  MTI  followed  by  incoherent 
integration. 


BACKGROUND 

This  report  describes  the  results  of  work  that  is  in  progress 
at  the  Naval  Ocean  Systems  Center  (NOSC),  San  Diego,  California,  and 
at  ORINCON  Corporation,  La  Jolla,  California.  Further  results  will 
be  described  in  an  expanded  version  of  this  report  which  will  be  sub- 
mitted for  publication  in  the  IEEE  Transactions  on^  Aerospace  and 
Electronic  Systems. 

INTRODUCTION 

The  output  sequence  from  a nonrecursive  MTI  is  a linear  combination 
of  the  inputs  from  successive  pulse  repetition  intervals  (PRI's).  Thus, 
the  sequence  at  the  MTI  output  is  correlated  even  though  the  input  se- 
quence may  be  uncorrelated.  When  detection  is  performed  by  incoherent 
integration  of  successive  MTI  outputs  (for  each  range  bin),  the  variables 
integrated  are  not  statistically  independent.  This  poses  two  problems. 
First,  the  determination  of  the  detection  threshold  must  take  into 
account  the  noise  correlation,  and  thus  graphs  or  tables  of  thresholds 
(e.g.,  Pachares'  table  [1])  previously  published  cannot  be  used.  Sec- 
ond, the  noise  correlation  degrades  detection  performance,  as  has  been 
indicated  in  the  literature  [2,3,4].  Trunk  [2]  considered  both  problems 
and  obtained  detection  thresholds  and  detection  sensitivity,  by  Monte  Carlo 
simulation,  for  a linear  integrator  at  the  output  of  a binary-weighted  MTI. 
Hall  and  Ward  [3]  and  Kretschmer  [4]  estimated  the  detection  sensitivity 
degradation  for  a square-law  integrator  by  computing  the  "effective  number" 
N of  pulses  from  a ratio  of  variances  and  squared  means. 


This  report  serves  two  purposes.  The  first  is  to  provide  an  exact 


computational  method  (and  the  associated  theory)  for  determining  the 
detection  threshold  to  be  used  when  a square-law  integrator  is  applied 
to  the  output  from  a nonrecursive  MTI.  Using  this  method,  a table  of 
bias  levels  is  computed  for  a two-pulse  (single-delay)  MTI  for  several 
values  of  the  number  N of  pulses  integrated  and  the  probability  PFA  of 
a false  alarm. 

The  second  purpose  is  to  derive  the  characteristic  function  of  the 
output  of  the  square-law  integrator  when  the  MTI  input  contains  a sig- 
nal with  Doppler  frequency  f^  and  amplitude  G.  While  inversion  of  this 
characteristic  function  (to  find  the  associated  probability  density) 
has  proven  to  be  intractable,  the  techniques  of  Helstrom  [5]  can  be  used 
to  estimate  the  detection  probability.* 

The  general  theory  is  given  for  obtaining  the  characteristic  func- 
tion of  the  sum  of  squared  magnitudes  of  correlated,  circular- 
Gaussian  [6]  random  variables.  (The  particular  circular-Gaussian  process 
considered  here  corresponds  to  the  in-phase  and  quadrature  components  of 
the  complex  envelope  of  a Gaussian  stochastic  process.)  For  the  case 
where  the  variables  have  zero  mean  (i.e.,  the  no-signal  case),  the 
characteristic  function  is  inverted,  and  the  probability  distribution  is 
obtained.  The  application  of  the  general  theory  to  the  computation  of  bias 


Estimates  of  the  detection  probability  will  be  given  in  the  journal 
publication  which  will  follow  at  a later  date.  Comparisons  will  then 
be  made  with  the  estimates  of  Hall  and  Ward  [3]  and  the  simulation 
resul ts  of  Trunk  [2] . 
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levels  for  a nonrecursive  MTI  is  discussed  in  detail.  Also,  applica- 
tions to  signal  detection  based  on  incoherent  integration  of  overlapping 
discrete  Fourier  transforms  (DFT's)  are  considered  briefly. 

The  theory  given  here  is  similar  to  derivations  by  other  authors, 
such  as  Kac  and  Siegert  [7],  Emerson  [8],  and  Meyer  and  Middleton  [9]. 

For  example,  Kac  and  Siegert  [7]  obtained  the  characteristic  function 
of  the  (unsampled)  output  from  a video  amplifier  which  follows  a 
square-law  detector,  and  expressed  their  results  in  terms  of  eigenvalues 
and  eigenfunctions  of  an  integral  equation.  Emerson  [8]  considered  the 
same  problem  and  outlined  procedures  for  determining  probability  density 
functions  directly  without  solving  the  eigenvalue  problem  or  inverting 
the  characteristic  function.  Meyer  and  Middleton  [9]  extended  the  re- 
sults of  Kac  and  Siegert  and  obtained  an  explicit  solution  for  the 
integral  equation  involving  the  autocorrelation  function  of  the  noise. 

Our  work  differs  from  these  previous  analyses  in  that  we  approach  the 
problem  from  a sampled-data  viewpoint  and  develop  an  exact  computational 
method  which  provides  a straightforward  means  of  analyzing  the  performance 
of  digital  detection  processors  which  employ  incoherent  integration  of 
correlated  random  variables. 

NONRECURSIVE  MTI 

For  each  range  bin,  the  output  from  a nonrecursive  MTI  is  a linear 
combination  of  the  inputs  from  L successive  pulse  repetition  intervals 
(PRI's).  That  is,  the  complex  MTI  output  E^  is* 

*The  output  is  a complex  variable  obtained  from  the  in-phase  and  quadrature 
components.  The  index  k denotes  the  PRI;  the  index  that  would  denote  the 
range  bin  has  been  suppressed. 
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(1) 


L 

Ek  = I “t  <V»L  = Xk  * ,Vk 
1=1 

where  W^,  l^,  • ••>  are  fixed  weights  which,  for  this  analysis,  are 
assumed  to  be  real.  The  real  and  imaginary  parts  of  the  input  sequence 
Gl*  G2’  •••»  are  assumed  to  be  independent  normally  distributed  random 
variables  with  variance  a2,  and  with  means  rrij,  m 2,  ...,  and  Cj,  C 2,  .... 
respectively.  Thus,  the  output  Ek  4xk  + iYk,  k=l,2,...  , is  such  that 
anc*  \ are  normally  distributed  with  means 

L 

uk  = E[Xk]  = 2 W£  mk-*+L  (2) 

1=1 


L 

Yk  = E^Yk^  zL  W£  Ck-£+L 
1=1 


THEORY 


With  Xj,  X^,  ....  X^  and  Yj,  •••>  Y^  defined  as  above,  their 
joint  probability  density  is  given  by 

f(x,z)  = (2tt)"^ | A|exp  |-  \ t(x-w)A(x,V)+(^-1)A(y.,-x,)l)  (4) 

where  A is  the  inverse  of  the  covariance  matrix  R determined  by  (3a); 

_x,  y.,  etc.,  denote  N-dimensional  row  vectors;  x.' > y.1 , etc.,  are  their 
transposes  (column  vectors);  and  |A|  is  the  determinant  of  the  matrix  A. 

From  Kendall  and  Stuart  [10,  page  347],  there  is  an  orthogonal 
matrix  (i.e.,  transformation)  T such  that  T'AT  = D,  where  D is  a 
diagonal  matrix.  Also,  the  elements  on  the  diagonal  of  D are  the  eigen- 
values of  A [11,  page  186].  Applying  this  transformation  T to  x and  y_ 
(by  letting  x'  = Tz/  and  y/  = Tw'),  the  new  variables  z and  w have  joint 
density 

f(z,w)  = (2tt)_N  | D|exp  {-|[(z-a)D(z,-al  ) + (w-b)D(w'-6i  )])  , (5) 


where  p_‘  = Ta'  and  x'  = Tjj' . 

Note  that  since  T is  orthogonal, 

N N 

S = ^ + yk  = +y-y-'  = IT'Tl'  + wT'Tw'  = 7 wk  + zk  - 

k=l  k=l 

(6) 
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From  (6)  and  the  discussion  in  the  previous  section,  it  follows  that  S 
can  be  interpreted  as  the  result  of  square-law  integration  of  N outputs 
from  a nonrecursive  MTI.  We  now  determine  the  characteristic  function 
of  S from  the  joint  density  (5). 

The  change  of  variables  zk  = rk  cos  ek  and  wk  = rk  sin  @k  in  (5) 
and  integration  with  respect  to  e^,  @2.  ....  0^  gives  the  joint  density 
of  rr  r2,  . , rN: 


f(r)  • n dkrk  exp  (rjj  * S*)'  WkS|c> 


2 A 2 2 

where  gk  = + ek>  the  dk's  are  the  eigenvalues  of  A,  and  I Q ( • ) is  the 

modified  Bessel  function  of  the  first  kind,  order  zero.  Finally,  letting 


qk  = rk,  the  joint  density  of  q ^ , q2> 


•••»  hn 


Pm  is 


dk  dk  / 2 

f (a)  = IT  exp  - -?r  (qk  + \ ) Io(9kdk  \ ) ■ 

k=l 


From  (8)  it  follows  that  q^,  q2,  . ..,  q^  are  independent  random  variables; 
also,  the  transformations  above  are  such  that 


s = I dk  ■ 


Thus,  the  characteristic  function  (ch.  f . ) of  S is. given  by 
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4>su)  - n 

k=l 

where  4^(10)  is  the  ch.  f.  of  q^.  From  Weinstock  [12,  page  177],  the 
ch.  f.  <t>k(w)  is 


;W=(1_^)-1  .«p  [,**/(  1 -ila) 


and,  therefore, 


n (i-x2)'1  exp  1 W/  (l-^) 

k=l  k k 


Inversion  of  (co ) provides  the  probability  density  f$(s)  of  the 

random  variable  S.  From  this  density  we  can  determine  the  false-alarm 

and  detection  probabilities  associated  with  a detection  system  which  is 

based  on  the  square-law  integration  of  the  outputs  from  a nonrecursive 

o 

MTI.  Inversion  of  (11)  for  the  signal -present  case  (i.e.,  when  g^  is 
nonzero)  has  proven  to  be  intractable.  However,  for  the  no-signal  case 
(g^  = 0)>  inversion  can  be  accomplished  as  described  in  the  next  sectior 


BIAS  LEVELS  FOR  SQUARE-LAW  INTEGRATION 

In  this  section  we  obtain  an  exact  computational  method  for  deter- 
mining the  bias  levels  (detection  thresholds)  to  be  used  when  a square- 
law  integrator  is  applied  to  the  output  from  a nonrecursive  MTI. 


t; 

t. 


When  no  signal  is  present  at  the  MTI  input,  the  means  of  and  Y| 
(the  real  and  imaginary  parts  of  the  MTI  output)  are  zero.  Thus,  from 
(11) , the  ch.  f . of  S is 


When  the  eigenvalues  d^,  d^,  . ..,  d^  are  all  distinct  (which  is  the  case 
for  the  covariance  matrices  considered  here),  (12)  can  be  inverted  to 
give 


fS(s ) 


* I 


dk  V 


-sdk/2 


k=l 


(13a) 


where 


Therefore,  the  false-alarm  probability  PFA  is  given  by 


PFA  = 


N 

I 

k=l 


Pke 


-Tdk/2 


(13b) 


(14) 


where  T is  the  detection  threshold.  From  (14)  it  follows  that  PFA  is 
uniquely  determined  by  the  threshold  T and  the  eigenvalues  of  A (which 
is  the  inverse  of  the  covariance  matrix  R).  However,  R is  the  matrix 
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that  is  given  (e.g.,  from  equation  (3))  and  its  inversion  can  be 
computationally  difficult.  Since  the  eigenvalues  of  A are  reciprocals 
of  the  eigenvalues  of  R,  it  is  not  necessary  to  invert  R.  That  is,  it 
is  only  necessary  to  find  the  eigenvalues  X^,  x^,  ....  x^  of  R.  (Note 
that  x^  = d^.)  1°  terms  of  the  eigenvalues  of  R,  PFA  becomes 


PFA  = 


N 

y 

k=l 


Bke 


"T/2x, 


where 


B 


k 


N 

n 

t=i 


(15a) 


(15b) 


The  computational  method  for  obtaining  the  bias  level  T^  is  as 
follows:  For  each  N,  the  eigenvalues  x^,  x^,  ....  x^  are  determined  and 
the  coefficients  B^,  B2,  B^  are  obtained  from  (15b).  For  each 
specified  value  of  PFA,  (15a)  is  solved  for  the  bias  level  T = T^  ppA 
using  Newton's  method  (or  another  numerical  technique).  The  computations 
are  illustrated  in  the  next  section  where  a particular  example  is 
discussed. 
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TWO-PULSE  MTI  FOLLOWED  BY  INCOHERENT  INTEGRATION 


For  a two-pulse  MTI  (L  = 2 in  (3))  the  covariance  matrix  R^ 
has  the  form 


2 2 2 2 

where  a = (W1  + W^)  a and  b = W1  W2  a . (The  subscript  N denotes  the 
dimension  of  R and  is  also  the  number  of  MTI  outputs  incoherently  inte- 
grated.) If  binomial  weights  [13]  are  used  (i.e.,  Wj^  = 1,  Wg  = -1), 

2 2 

then  a = 2a  and  b = -a  . In  the  following,  a two-pulse  MTI  with  binomial 

weights  is  assumed.  Also,  since  the  detection  thresholds  are 

2 

normalized  by  the  noise  variance  a , no  loss  of  generality  occurs  if 
2 

we  assume  o =1. 

Bias  levels 

To  obtain  the  eigenvalues  of  R^,  we  compute  the  determinant  of  the 
matrix  R^  - Xl^  and  find  the  roots  of  the  polynomial  equation 

|Rn  - AIN|  = 0.  (17) 


For  the  case  being  considered,  the  polynomials  satisfy  a recursive 
relation  as  a function  of  N.  We  make  the  change  of  variable  u = 2-x  in 
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(17)  and  denote  by  DN(u)  the  polynomial  in  u obtained  from  (17).  The 
following  recursive  formula  is  easily  verified: 


Dj(u)  = u 
D2(u)  = u2  - 1 

(18) 

dn(u)  = u Dn1(u)  - Dn_2(u) 

N = 3,4,... 

From  (18)  the  derivative  D^(u)  is  also  obtained  recursively: 

D‘(u)  = 1 
D£(u)  = 2u 

(19) 

DN(u)  = u DN-1(u)  " DN-2(u)  + DN-1^u) 

N = 3,4,... 

Using  (18)  and  (19),  Newton's  method  can  be  applied  to  compute  the 
roots  of  (17). 


Equations  (18)  and  (19)  were  derived  under  the  assumption  that 
a = 2 and  b = -l  in  (16).  However,  it  is  easily  shown  that  the  solutions 


of  (17)  in  terms  of  the  variable  u defined  above  provide  the  most  gen- 
eral solutions  for  any  matrix  of  the  form  defined  in  (16).  In  (17),  if 
we  let  v=(a-A)/|b|,  then  the  resulting  polynomial  in  v is  proportional 
to  D^U),  where  0^(0  is  defined  in  (18).  Thus,  the  roots  u^, 
k = 1,2,  ...  , N,  of  (18)  can  be  used  to  obtain  the  eigenvalues  of  any 
matrix  of  the  form  (16)  by  defining  A^  = a - | b ( u ^ . 

To  illustrate  the  location  of  the  eigenvalues  of  RN>  Figure  1 is  a 
graph  of  Aj^,  ....  A^  as  a function  of  N for  N between  1 and  20.  The 
continuous  curves  are  obtained  by  connecting  the  mth  largest  eigenvalue 
for  each  N.  Because  of  symmetry  about  A = 2,  only  the  eigenvalues  in  the 
interval  [2,4]  are  shown.  (The  eigenvalues  for  a particular  value  of  N 
are  found  in  Figure  1 at  the  intersection  of  the  corresponding  horizon- 
tal line  with  the  continuous  curves.) 

Table  1 gives  the  bias  levels  required  to  obtain  the  indicated 
values  of  PFA  for  values  of  N between  2 and  39.  Equations  (18)  and  (19) 
were  used  to  find  the  eigenvalues  and  (15)  was  used  to  determine  T by  the 
method  described  in  the  previous  section. 

Table  1 was  computed  for  the  case  where  a = 2 and  b =-l  in  (16).  How- 
ever, the  results  are  more  generally  applicable.  From  the  discussion 
above,  the  solutions  of  (18)  can  be  used  to  determine  the  eigenvalues  A^ 
from  Ak  = a - |b|uk.  If  a = 2|b|,  then  A^  = |b|(2-uk)  and  the  required 

threshold  is  |b|T,  where  T is  the  threshold  given  in  Table  1.  For  ex- 
2 

ample,  if  a t 1 (as  was  assumed  previously)  the  required  threshold  is 
2 2 

Ta  , where  a is  the  variance  of  the  real  and  imaginary  inputs  to  the  MT I 
It  should  be  noted  that  comparisons  of  the  bias  levels  in  Table  1 
with  Pachares’  table  [1]  must  take  into  account  a difference  in  normalization 
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1 

n 

3 

4 

5 

6 

n 

16.2394 

30.0637 

43.8793 

57.6948 

71.5103 

85.3258 

3 

22.5577 

38.6635 

54.4551 

70.1904 

85.9148 

101 .6369 

4 

28.4521 

46.2636 

63.3186 

80.1319 

96.8537 

113.5390 

5 

34.0875 

53.3038 

71.3648 

88.9887 

106.4079 

123.7218 

6 

39.5437 

59.9806 

78.9033 

97.2102 

115.2028 

133.0168 

7 

44.8678 

66.3986 

86.0883 

104.9981 

123.4918 

141 .7368 

8 

50.0896 

72.6211 

93.0118 

112.4719 

131.4218 

150.0577 

9 

55.2293 

78.6892 

99.7301 

119.7009 

139.0738 

158.0717 

10 

60.3013 

84.6312 

106.2821 

126.7323 

146.5028 

165.8410 

11 

65.3162 

90.4679 

112.6956 

133.5999 

153.7474 

173.4087 

12 

70.2S22 

96.1920 

118.9632 

140.2953 

160.7979 

180.7633 

13 

75.2057 

101.8642 

125.1612 

146.9088 

167.7576 

188.0205 

14 

80.0919 

107.4673 

131.2694 

153.4165 

174.5986 

195.1484 

15 

84.9448 

113.0105 

137.2983 

159.8306 

181.3347 

202.1620 

16 

89.7680 

118.5126 

143.2727 

166.1804 

187.9989 

209.0978 

17 

94.5644 

123.9505 

149.1648 

172.4335 

194.5545 

215.9149 

18 

99.3363 

129.3446 

154.9998 

178.6193 

201 .0347 

222.6498 

19 

104.0860 

134.6885 

160.7703 

184.7295 

207.4302 

229 . 2922 

20 

108.8152 

140.0062 

166.5055 

190.7978 

213.7786 

235.8837 

21 

113.5255 

145.2897 

172.1961 

196.8136 

220.0682 

242.4109 

118.2183 

150.5416 

177.8454 

202.7807 

226.3031 

248.8786 

23 

122.8948 

155.7640 

183.4564 

208.7026 

232.4873 

255.2910 

24 

127.5562 

160.9591 

189.0316 

214.5822 

238.6241 

261.6516 

25 

132.2034 

166.1286 

194.5734 

220.4224 

244.7165 

267.9638 

26 

136.8373 

171.2741 

200.0837 

226.2256 

250.7673 

274.2305 

27 

141.4586 

176.3969 

205.5646 

231.9939 

256.7788 

280.4545 

28 

146.0681 

181.4984 

211.0176 

237.7292 

262.7534 

286.6379 

29 

150.6665 

186.5798 

216.4443 

243.4335 

268.6930 

292.7831 

30 

155.2543 

191.6422 

221.8461 

249.1084 

274.5995 

298.8921 

31 

159.8320 

196.6907 

227.2291 

254.7607 

280.4805 

304.9732 

32 

164.4002 

201.7135 

232.5798 

260.3755 

286.3194 

311.0083 

33 

163.9593 

206.7242 

237.9140 

265.9704 

292.1358 

317.0188 

34 

173.5098 

211.7193 

243.2279 

271.5412 

297.9251 

322.9994 

35 

178.0520 

216.6995 

248.5222 

277.0889 

303.6882 

328.9515 

36 

182.5862 

221.6654 

253.7978 

282.6145 

309.4265 

334.8764 

37 

187.1129 

226.6177 

259.0556 

288.1189 

315.1409 

340.7750 

38 

191.6324 

231.5570 

264.2963 

293.6031 

320.8324 

346.6485 

39 

196.1449 

236.4837 

269.5205 

299.0677 

326.5019 

352.4979 

Table  1.  Bias  levels  for  square-law  integration  of  N outputs 
from  a binary-weighted,  two-pulse  nonrecursive  MTI. 
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7 

8 

9 

10 

11 

12 

n 

99.1414 

112.9569 

126.7724 

140.5879 

154.4034 

168.2189 

3 

117.3586 

133.0801 

148.8014 

164.5227 

180.2440 

195.9652 

4 

130.2093 

146.8733 

163.5347 

180.1950 

196.8548 

213.5144 

5 

140.9789 

158.2045 

175.4125 

192.6106 

209.8029 

226.9920 

6 

150.7227 

168.3606 

185.9548 

203.5204 

221.0672 

238.6014 

7 

159.8229 

177.8028 

195.7096 

213.5649 

231.3835 

249.1755 

8 

168.4865 

186.7724 

204.9561 

223.0651 

241.1184 

259. 1296 

9 

176.8173 

195.3849 

213.8226 

232.1629 

250.4288 

268.6371 

10 

184.8843 

203.7164 

222.3920 

240.9486 

259.4128 

277.8040 

11 

192.7346 

211.8180 

230.7197 

249.4816 

268.1337 

286.6982 

12 

200.3552 

219.6747 

238.7885 

257.7427 

276.5705 

295.2966 

13 

207.8735 

227.4256 

246.7489 

265.8937 

284.8961 

303.7831 

14 

215.2534 

235.0303 

254.5564 

273.8856 

293.0571 

312.1002 

15 

222.5113 

242.5062 

262.2291 

281.7376 

301.0736 

320.2684 

16 

229.6866 

249.8958 

269.8125 

289.4980 

308.9967 

328.3418 

17 

236.7341 

257.1496 

277.2529 

297.1086 

316.7637 

336.2535 

18 

243.6940 

264.3109 

284.5967 

304.6191 

324.4275 

344.0591 

19 

250.5545 

271.3669 

291.8298 

312.0140 

331.9713 

351.740 7 

20 

257.3611 

278.3665 

299.0047 

319.3492 

339.4543 

359.3606 

21 

264.0989 

285.2936 

306.1036 

326.6057 

346.8559 

366.8970 

nn 

270.7730 

292.1534 

313.1322 

333.7890 

354.1821 

374.3557 

23 

277.3878 

298.9505 

320.0952 

340.9042 

361.4378 

381.7420 

24 

283.9473 

305.6891 

326.9969 

347.9557 

368.6278 

389.0607 

25 

290.4549 

312.3729 

333.8412 

354.9475 

375.7561 

396.3159 

26 

296.9138 

319.0051 

340.6315 

361 .8832 

382.8263 

403.5113 

27 

303.3267 

325.5887 

347.3709 

368. 7659 

389.8418 

410.6504 

28 

309.6963 

332.1264 

354.0622 

375.5985 

396.8054 

417.7361 

29 

316.0248 

338.6207 

360.7079 

382.3837 

403.7199 

424.7712 

30 

322.3144 

345.0737 

367.3103 

389.1238 

410.5878 

431.7582 

31 

328.5741 

351 .4951 

373.8797 

395.8297 

417.4204 

438.7090 

32 

334.7844 

357.8640 

380.3937 

402.4776 

424.1926 

445.5973 

33 

340.9682 

364.2049 

386.8784 

409.0951 

430.9334 

452.4534 

34 

347.1200 

370.5118 

393.3274 

415.6753 

437.6357 

459.2697 

35 

353.2411 

376.7862 

399.7422 

422.2199 

444.3010 

466.0478 

36 

359.3328 

383.0294 

406.1243 

428.7303 

450.9308 

472.7893 

37 

365.3964 

389.2427 

412.4750 

435.2080 

457.5267 

479.4957 

38 

371.4329 

395.4273 

418.7954 

441 .6541 

464.0898 

486.1684 

39 

377.4435 

401.5843 

425.0868 

448.0699 

470.6215 

492.8086 

Table  1.  (Continued) 
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To  obtain  the  ch.  f.  <f><.(w)  from  (11),  we  must  determine  + 

where  ak  and  Bk  are  components  of  the  vectors  obtained  from  the  trans- 
formations jV  = Ta‘  and  x'  = T$ 1 . It  can  be  shown  [11]  that  the 
matrix  T has  rows  which  are  the  unit-length  eigenvectors  ^ corres- 
ponding to  the  eigenvalues  Ak>  k=  1,2,  ...»  N.  Therefore,  is  given 
by 

9k2  = (t*  u')2  + x')2  (21) 

The  eigenvectors  t^  are  obtained  from  the  relationship 

RN  = Xk  (22) 

and,  for  the  covariance  matrix  given  by  (16)  (with  a = 2 and  b=-l), 
equation  (22)  leads  to  the  relationship 

tkj  = Dj-1^2  " Ak)  tkl  (23a) 

j = 2,3,  ...  , N 

where 

h - ^kl’  tk2’  ’ W (23b) 

and  D j ( • ) is  defined  above  in  (17).  In  (23a),  tkl  must  be  chosen  such 
that  t^  has  unit  length.  Thus,  we  can  define  the  eigenvector  t^  from 
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2 2? 

To  obtain  the  ch.  f.  <t»^(to)  from  (11),  we  must  determine  4 + 6k  ’ 

where  and  ek  are  components  of  the  vectors  obtained  from  the  trans- 
formations p_'  = Ta'  and  x'  = Te/ . It  can  be  shown  [11]  that  the 

matrix  T has  rows  which  are  the  unit-length  eigenvectors  corres- 

2 

ponding  to  the  eigenvalues  k=l,2,  ....  N.  Therefore,  gk  is  given 
by 

\ = (t^  u')2  + (t^  x')2  (21) 

The  eigenvectors  t^  are  obtained  from  the  relationship 

RN  H.  = Xk  ^k  (22> 


and,  for  the  covariance  matrix  given  by  (16)  (with  a = 2 and  b=-l), 
equation  (22)  leads  to  the  relationship 

‘kj  ■ Dj-i(2-»k)  ‘u  <23a> 

j ■ 2,3 N 

where 

h - ^kl’  tk2 ’ ’ W (23b) 

and  Dj(-)  is  defined  above  in  (17).  In  (23a),  tkj  must  be  chosen  such 
that  has  unit  length.  Thus,  we  can  define  the  eigenvector  t^  from 


tkl  = Ek 


lkj " Ek ' Dj-1(2-V 


j = 2,3,  ...  , N 


(24a) 


where 


' N-l 

1 + y D2(2-aJ 


(24b) 


Now,  rewriting  (21),  we  have 


2 _ 2 , 2 
gk  ' ak  + 6k 


N N 


X I K + >m)  hi  h 


(25a) 


£=1  m=l 


4G2  sin2  (tt  fd/PRF) 


N 

y t, 


(25b) 


with  (25b)  obtained  from  (25a)  using  the  relationship  defined  in  (20). 
Equation  (25b)  can  be  further  simplified.  Note  that 


I ‘kt  = i ik  ■ 


where  i 4 (1,1,  ....  l),  and  from  (22): 


i»N  - \ i y ■ »k  X 
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But  i_  = (1,0,  ....  0,1),  and  (26)  reduces  to 


^ tkd  ^kl  + W/Xk 


n=l 


Using  (27),  (25b),  and  (24)  we  have,  finally. 


(27) 


2 4G2  sin2  (tt  fd/PRF)  / 1 + D^  (2  - Afc)\ 2 


N-l 

1 * X Df(2  - V 

£=1 


(28) 


APPLICATIONS  TO  OVERLAPPING  DFT'S 

The  discussion  above  pertained  to  the  square-law  integration  of  the 
outputs  from  a nonrecursive  MTI.  We  now  consider  briefly  the  applica- 
tion of  the  theory  to  square-law  integration  of  coefficients  of  over- 
lapping DFT's. 

Signal  detection  is  often  accomplished  by  computing  the  discrete 
Fourier  transform  (DFT)  of  a sequence  GQ,  G^,  . . . , G^  ^ of  sensor  out- 
puts and  comparing  the  magni tudes-squared  of  the  DFT  coefficients  with 
a threshold.  The  length  M of  the  DFT  is  determined  by  several  factors, 
including  the  required  signal  resolution  in  the  frequency  domain  [14]. 

To  improve  detection,  it  is  sometimes  necessary  to  perform  incoher- 
ent integration  of  the  DFT  outputs.  That  is,  a sequence  of  DFT's  is 
computed  and  the  magni tudes-squared  of  successive  coefficients  are 
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summed  and  compared  with  a threshold.  The  sequence  of  DFT's  can  be  con- 
tiguous, as  indicated  in  Figure  3a,  or  overlapping,  as  indicated  in  Fig- 
ure 3b.  When  the  DFT's  are  nonoverlapping,  the  quantities  incoherently 
integrated  are  independent;  however,  when  the  DFT's  are  overlapping,  the 
integration  involves  dependent  (correlated)  random  variables.  (Since  we 
are  dealing  exclusively  with  Gaussian  processes,  correlation  and  depend- 
ence are  used  interchangeably.) 

The  correlation  between  coefficients  of  successive  transforms  de- 
pends both  on  the  degree  of  overlap  and  on  the  "window-function"  used  in 
the  DFT  [14,  page  56].  The  qth  coefficient  of  the  pth  DFT  is  defined  as 


M-l 

^qp  ^k  ^k+p(M-r) 

k=0 


e-2niqk/M 


p = 0, 1,  . . . ; 0<r<M-l 


(29) 


N-l 

where  { W.  } is  the  window  function,  and  the  parameter  r defines  the 
K k=0 

degree  of  overlap.  (Note  that  r is  defined  differently  here  than  by 
Harris  [14,  page  56].)  Let  the  real  and  imaginary  parts  of  Fqp  be 
defined  by 


+ i Y 


qp 


(30) 


and  assume  that  the  inputs  G^,  k = 0,l,  ...  , have  real  and  imaginary 

2 

parts  that  are  1 1 D Gaussian  random  variables  with  variance  a . Then, 
the  covariance  of  Xqp,  Xq>p+u  (and  Yqp,  Yq>p+u)  is  given  by 
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4 

j 


i 


DFT  1 DPT  2 DPT  3 


Figure  3a.  Partitioning  of  input  sequence  for  nonoverlapping  DFT's. 


—►I  M-r 


DPT  1 


DPT  2 


DFT  3 


Figure  3b.  Partitioning  of  input  sequence  for  overlapping  DFT’s. 


* 


M- 1+u ( r-M) 

*u  - »2  I Vu(M-r)  • 0 i u (31) 

n=o 

= 0 , otherwise. 

(The  index  q has  been  suppressed  in  (31)  since  the  covariance  is  inde- 
pendent of  q.)  Using  (31),  the  covariance  matrix  of  the  real  and  imag- 
inary parts  of  Fqg,  F^,  ...  , Fq  ^ can  be  determined  and  the  previ- 
ously derived  theory  can  be  applied  to  determine  the  bias  levels  to  be 
used.  If  the  presence  of  a signal  is  represented  by  having  nonzero- 
mean  inputs  (such  as  defined  by  equation  (20)  in  the  previous  section), 
the  characteristic  function  for  the  signal -present  case  can  also  be 
derived. 

If  the  signal  to  be  detected  is  a narrowband  Gaussian  process,  the 
presence  of  such  a signal  at  the  OFT  input  results  in  an  increase  in  the 
variance  of  the  real  and  imaginary  parts  of  the  DFT  coefficient  which 
corresponds  to  the  center  frequency  of  the  input  process.*  Thus,  for 
this  coefficient  the  covariance  when  a signal  is  present  is 


♦Because  of  spectral  leakage  [14,  page  52],  the  signal  will  increase  the 
variance  of  all  DFT  coefficients.  We  assume  that  the  bandwidth  of  the 
input  process  is  narrow  enough  that  only  the  DFT  coefficient  correspond- 
ing to  the  center  frequency  of  the  process  will  have  a significant 
response. 


1; 


23 


2 

where  o<-  is  the  variance  of  the  signal  component  of  the  DFT  input,  and 

4>u  is  the  noise-only  covariance  given  by  (31).  If  the  SNR  t| > is  defined 

2 2 I 

such  that  <i i = /o  , the  covariance  $u  <.  is 

^u,S  = (1  + *>  *u  ' (33)  . I 

It  follows  that  the  probability  of  detection  Pd  can  be  computed  from 

(15a)  with  T replaced  by  T/ ( 1 + ip ) . 

As  an  example,  for  50%  overlap  the  covariance  matrix  has  the  form 

of  (16)  with  a = $o  and  b = (j> ^ , and  recursive  formulas  similar  to  (18)  and 

(19)  can  be  derived.  Fora  rectangular  window  (i.e.,  Wk=l,  k = 0,l,  ...  , 

2 M2 

M-  1),  a = Mo  and  b = ^ o . Thus,  the  thresholds  tabulated  in  Table  1 are 

2 

applicable  if  T is  replaced  by  MTo  /2.  This  follows  from  the  discussion 
of  Table  1 above  and  the  fact  that  a = 2b. 

CONCLUSIONS  AND  DISCUSSION 

The  general  theory  has  been  given  for  obtaining  the  ch.  f.  of  the 
sum  of  magnitudes  squared  of  correlated  circular-Gaussian  random  variables. 

This  theory  has  been  applied  to  the  case  of  incoherent  integration  of  the 
outputs  from  a nonrecursive  MTI.  In  particular,  the  ch.  f.  has  been  de- 
termined explicitly  for  a two-pulse  MTI,  and  the  threshold  levels  which 
provide  specified  false-alarm  probabilities  have  been  determined  and 
tabulated  for  several  values  of  the  number  N of  MTI  outputs  integrated. 

The  application  of  the  theory  to  incoherent  integration  of  the  co- 
efficients of  overlapping  DFT's  has  also  been  shown.  It  has  been  shown  that 


' 


t 


when  the  DFT's  overlap  by  50%,  the  tabulated  thresholds  for  the  two- 
pulse  MTI  can  be  used  after  appropriate  scaling. 


Detection  probabilities  have  not  been  computed  since  inversion  of 
the  ch.  f.  for  the  signal -present  case  has,  thus  far,  proven  to  be  intract- 
able. However,  the  technique  of  Hel strom  [5]  will  be  applied  to  estimate 
the  detection  probability  and  the  results  will  be  published  later.  These 
results  will  be  compared  with  the  estimates  of  Hall  and  Ward  [3]  and 
Trunk  [2],  and  the  degradation  in  performance  due  to  noise  correlation 
will  be  evaluated. 
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